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The accelerating expansion of our universe at present could be driven by an unknown energy 
component (Dark Energy) or a modification of general relativity (Modified Gravity). In this note 
we revisit the constraints on a phenomenological model which interpolates between the pure ACDM 
model and the Dvali-Gabadadze-Porrati (DGP) braneworld model with an additional parameter a. 
Combining the cosmic microwave background (CMB), baryon acoustic oscillations (BAO) and type 
la supernovae (SNIa), as well as some high-redshift observations, such as the gamma-ray bursts 
(GRB) and the measurements of linear growth factors (LGF), we obtain the tight constraint on the 
parameter a = 0.254 ± 0.153 (68% C.L.), which implies that the flat DGP model is incompatible 
with the current observations, while the pure ACDM model still fits the data very well. Finally, we 
simulate the future measurements with higher precisions and find that the constraint on a can be 
improved by a factor two, when compared to the present constraints. 
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I. INTRODUCTION 

Current cosmological observations, such as the CMB 
measurements of temperature anisotropies and polariza- 
tion at high redshift z ~ 1090 [l[ and the redshift- 
distance measurements of SNIa at z < 2 , have demon- 
strated that the universe is now undergoing an acceler- 
ated phase of expansion. The simplest explanation is 
that this behavior is driven by the cosmological constant 
or the dynamical dark energy models, which suffers from 
the severe coincidence and fine-tuning problems [![. On 
the other hand, this observed late-time acceleration of 
the expansion on the large scales could also be caused 
by some modifications of general relativity, such as the 
scalar-tensor 0] and f(R) theories Q, and gravitational 
slip @. 

One of the well-known examples is the DGP 
braneworld model 0, in which the gravity leaks off the 
four dimensional brane into the five dimensional space- 
time. On small scales gravity is bound to the four di- 
mensional brane and the general relativity is recovered 
to a good approximation. In the framework of flat DGP 
model, the Friedmann equation will be modified as [f| : 
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where r c — (Hq(1 — f2 m )) _1 is the crossover scale. At 
early times, Hr c 3> 1, the Friedmann equation of general 
relativity is recovered, while in the future, H — > = 
l/r c , the expansion is asymptotically de Sitter. Recently 
there have been a lot of interests in the phenomenological 
studies relevant to the DGP model in the literature [9j, 

In this note we investigate an interesting phenomeno- 
logical model, first introduced in Ref . [TlT| . which inter- 
polates between the pure ACDM model and the DGP 
model with an additional parameter a and presents the 



tight constraints from the current observations and fu- 
ture measurements. The paper is organized as follows: 
In Sec. II we describe the general formalism of the mod- 
ified gravity model. Sec. Ill contains the current obser- 
vations we use, and Sec. IV includes our main global 
fitting results. In Sec. V we present the forecasts from 
the future measurements, while Sec. VI is dedicated the 
conclusions. 



II. GENERAL FORMALISM 

In this phenomenological model, assuming the flatness 
of our universe, the Friedmann equation is modified as 
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where r c = i? _1 /(l — S7 m ) a ~ 2 . Thus, we can straight- 
forwardly rewrite the above equation and obtain the ex- 
pansion rate as following: 
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where the last term denotes the modification of the Fried- 
mann equation of general relativity: 
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Furthermore, we can obtain the effective equation of 
state: 

1 d\n5H 2 a, ,E'(z) , . 

^^- 1+ 3^n(lTI)^ 1+ 3 (1 + ^' (5) 

where the prime denotes the derivative with respect to 
the redshift z. 

In FigQ]we illustrate the evolutions of the effective en- 
ergy density £l a (z) = 1 — Ct m (z) and w e ff(z) for different 
values of parameter a. During the matter dominated era, 
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FIG. 1: The evolutions of effective energy density fl a (z) and 
effective equation of state w e s(z) for different values of a. 
The black solid lines are for a = (ACDM), the red dashed 
lines for a = 0.25 and the blue dotted lines for a — 1 (DGP). 
And the vertical line denotes today (z — 0). 



FIG. 2: The evolutions of linear growth g(a) = 8(a) /a for 
different values of a. The black solid lines are for a — 
(ACDM), the red dashed lines for a = 0.25 and the blue 
dotted lines for a = 1 (DGP). 



E(z) varies as (1 + z) 3 ^ 2 , which corresponds to the effec- 
tive equation of state: w c q = —1 + a/2. In the future 
z — > 0, with the matter density p m oc (1 + z) 3 — > 0, we 
have Weff(z) — > — 1 and ft a (z) — > 1. 

Besides the modification of the Friedmann equation, 
the flat DGP model also changes the growth function 
of density perturbation 8(a). Under assumptions of a 
quasi-static regime and sub-horizon scales, the correct 
evolution of perturbation was found [jj U3] : 

8+2H8- 4ttG (\ + p m 8 = , (6) 

where the dot denotes the derivative with respect to the 
cosmic time t, and the (3 factor is given by: 
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However, this phenomenological model Eq.© is a 
parametrization, so the situation is more complicated. 
One of the possible methods was found by Ref.[13j. In 
order to obtain the growth function of density perturba- 
tion within a covariant theory, the authors introduced a 
correction term and assumed the structure of modified 
theory of gravity to determine this term. Based on those 
assumptions, it was consequently found that the (3 factor 
was: 
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In the following analysis, we will use Eq.© and Eq.© 
to calculate the growth of density perturbation. 



Defined the normalized growth g(a) = 8(a) /a, the 
growth function Eq.([B]) can be rewritten as: 
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where the variable X(a) is the ratio of the matter density 
to the effective energy density X(a) = Q m (a)/Q a (a). In 
Figj2]we plot the linear growth factor g(a) as function of 
scale factor a for different values of a. One can see that 
the linear growth factor has been suppressed obviously as 
long as a is larger than zero. Thus, in the literature the 
linear growth has been widely used to study the modified 
gravity models, especially the DGP model 

Furthermore, the growth factor can be parameterized 
as [H: 
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where 7 is the growth index. And then the growth func- 
tion becomes: 
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For the pure ACDM model, the theoretical value of 7 is 
6/11 w 0.545, while 7 = 11/16 = 0.6875 in the flat DGP 
model [3. 

In the framework of this phenomenological model, we 
can easily see that the pure ACDM model and flat DGP 
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FIG. 3: The evolution of the (3 term (3/3) 1 when fixing a = 
-2. 



the universe: 0.1 < Sl m < 0.5. Furthermore, we con- 
strain the Hubble parameter to be uniformly in 4cr Hub- 
ble Space Telescope (HST) region: 0.4 < h < 1.0. The 
resulting plots are produced with CosmoloGUI 1 . 

In this section we will list the cosmological observations 
used in our calculations: CMB, BAO and SNIa measure- 
ments, as well as some high-redshift observations, such as 
the GRB and LGF data. We have taken the total likeli- 
hood to be the products of the separate likelihoods (£j) 
of these cosmological probes. In other words, defining 
X 2 L ,i = -2 log A, we get: 

2 _ 2 i 2 i 2 

AL,total — Xl,CMB "+" Xl.BAO "+" Xi,SNIa 

+ xi,GRB + Xl.LGF • ( 14 ) 

If the likelihood function is Gaussian, Xl coincides with 
the usual definition of \ 2 U P to an additive constant cor- 
responding to the logarithm of the normalization factor 
of £. 



model can be recovered when a — and a = 1, respec- 
tively. In order to be consistent with the cosmological 
observations, the a term should be very small in the early 
times, such as the Big Bang Nucleosynthesis (BBN) era 
(z ~ 10 9 ). This limit corresponds to the upper bound: 
a < 2 [1JJ. On the other hand, when a < 0, the ef- 
fective equation of state will become smaller than —1, 
which leads to the instability of linear growth of density 
perturbation due to the appearance of term Eq. © . 

At early times, such as the matter dominated era, we 
have E{z) oc (l + z) 3/2 and H/H 2 ~ -3/2. Thus, Eq.® 
becomes 
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since (1 + z) 3 > 1 at 2 < z < 1000. By contrast, at late 
times the matter energy density p m oc (1 + z) 3 — > and 



the expansion is asymptotically de Sitter, H 
then we have 
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Based on Eq. (fT2f and Eq. ([T3|) . we can straightforwardly 
see that as long as a < 0, during the evolution of universe 
the value of (3 should change the sign at one pivot redshift 
z t , which leads to j3\ Zt — and (3/3) 1 1 ^± ~~ * °°- I n 
Fig|3]we have shown the evolution of (3/3) _1 when fixing 
a = —2. There is an obvious singularity at a ~ 0.595. 
Therefore, based on these discussions above, we use a 
tophat prior on a as < a < 2 in our calculations. 



III. METHOD AND DATA 

In our calculations we assume a flat space and use a 
uniform prior on the present matter density fraction of 



A. CMB Data 

CMB measurement is sensitive to the distance to the 
decoupling epoch via the locations of peaks and troughs 
of the acoustic oscillations. Here we use the "WMAP 
distance information" obtained by the WMAP group [l| , 
which includes the "shift parameter" R, the "acoustic 
scale" I a, and the photon decoupling epoch z 2 . R and 
I a correspond to the ratio of angular diameter distance 
to the decoupling era over the Hubble horizon and the 
sound horizon at decoupling, respectively, given by: 



R = 



Ia = 



C 
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(15) 
(16) 



where x( z *) an d Xs{ z *) denote the comoving distance to 
z* and the comoving sound horizon at z„, respectively. 
The decoupling epoch z* is given by [l7l ]: 



z* = 1048[1 + 0.00124(lVi 2 )-°- 738 ][l + gi (n m h 2 ) 92 } 



where 
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(18) 



1 URL: http://www.sarahbridle.net/cosmologui/ 

2 In the revised version of WMAP5 paper [l|, they also extend the 
baryon density Q b h 2 into the WMAP distance information. But 
our calculations are not sensitive to Qi,h 2 and they also claim 
that this extension does not affect the constraints. Thus, we fix 
n b h 2 = 0.022765 to be the best fit value obtained by the WMAP 
group. 
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Table I. Inverse covariance matrix for the WMAP distance 
information I a, R and z*. The maximum likelihood values 
are R = 1.710, l A = 302.10 and z* = 1090.04, respectively. 
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We calculate the likelihood of the WMAP distance infor- 
mation as follows: 



1 >3 \ X j 



„dat£ 



), 



(19) 



where x = (R, Ia,z*) is the parameter vector and (C _1 )y 
is the inverse covariance matrix for the WMAP distance 
information shown in Tabic I. 



BAO Data 



The supernovae data we use in this paper are the recently 
released Union compilation (307 sample) from the Su- 
pernova Cosmology project [2], which include the recent 
samples of SNIa from the SNLS and ESSENCE survey, as 
well as some older data sets, and span the redshift range 
< z < 1.55. In the calculation of the likelihood from 
SNIa we have marginalized over the nuisance parameter, 
the absolute magnitude M, as done in Ref . [221] : 



where 
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(26) 



The BAO information has been already detected in the 
current galaxy redshift survey. The BAO can directly 
measure not only the angular diameter distance, D^(z), 
but also the expansion rate of the universe, H(z). But 
current BAO data are not accurate enough for extract- 
ing the information of Da(z) and H(z) separately [l8| . 
Therefore, one can only determine the following effective 
distance fl9j ]: 
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In this note we use the gaussian priors on the distance 
ratios r s (z d ) / D v (z): 



r s (z d )/D v {z = 0.20) 
r s (z d )/D v {z = 0.35) 



0.1980 ±0.0058 , 
0.1094 ± 0.0033 , (21) 



with a correlation coefficient of 0.39, extracted from the 
SDSS and 2dFGRS surveys where r s {z d ) is the co- 
moving sound horizon size and z d is the drag epoch at 
which baryons were released from photons given by [2lj : 
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where 



(22) 



(23) 



61 = 0.313(r! m /i 2 )- - 419 [l + 0.607(n m /i 2 ) - 674 ] 

6 2 = 0.238(n m h 2 ) 223 . 

C. SNIa Data 



The SNIa data give the luminosity distance as a func- 
tion of redshift 
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D. GRB Data 

GRBs can potentially be used to measure the luminos- 
ity distance out to higher redshift than SNIa. Recently, 
several empirical correlations between GRB observables 
were reported, and these findings have triggered intensive 
studies on the possibility of using GRBs as cosmological 
"standard" candles. However, due to the lack of low- 
redshift long GRBs data to calibrate these relations, in 
a cosmology-independent way, the parameters of the re- 
ported correlations are given, assuming an input cosmol- 
ogy, and obviously they depend on the same cosmological 
parameters that we would like to constrain. Thus, apply- 
ing such relations to constrain cosmological parameters 
leads to biased results. In Ref.[23[ the circular problem 
is naturally eliminated by marginalizing over the free pa- 
rameters involved in the correlations; in addition, some 
results show that these correlations do not change signif- 
icantly for a wide range of cosmological parameters [241 ] . 
Therefore, in this paper we use the 69 GRBs sample over 
a redshift range from z = 0.17 — 6.60 published in Ref.[25[ 
but we keep in mind the issues related to the "circular 
problem" that are more extensively discussed in Ref. [23[ . 



E. LGF Data 

As we point out above, the linear growth factor will 
be suppressed in the modified gravity model. It will be 
helpful using the measurements of linear growth factor 
to constrain the modified gravity models. Therefore, in 
Table II we list linear growth factors data we use in our 
analysis: the linear growth rate f = f2X, from galaxy 
power spectrum at low redshifts [26|, [2?], HH, H^, [3(| and 
lyman-cv growth factor measurement obtained with the 
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TABLE II. The currently available data for linear growth 



rates / 


we use 


in our 


analysis. 
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/ 


a 


Ref. 


0.15 


0.51 


0.11 


[26] 


0.35 


0.70 


0.18 


27] 


0.55 


0.75 


0.18 


[28] 


0.77 


0.91 


0.36 


[29] 


1.40 


0.90 


0.24 


[30] 


3.00 


1.46 


0.29 


[31] 



TABLE III. Constraints on the parameters a, Q m and 7. 
Here we have shown the mean values and errors from the 
current observations and the standard derivations from the 
future measurements. 





Q 
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CMB+BAO+SN 


0.263 ±0.175 


0.276 ±0.018 




All Real Data 


0.254 ±0.153 


0.277 ±0.017 


0.570 ±0.205 


Future 


0.07 


0.005 


0.050 



lyman-a power spectrum at z — 3 31]. It is worth rioting 
that the data points in Table II are obtained with assum- 
ing the ACDM model, thus, one should use these data 
very carefully, especially for the points obtained from 
Refs.p7l. [H, H3]. The corresponding \ 2 1S simply given 
by: 



* 2 = E 



^yth ydata^2 



(27) 



IV. NUMERICAL RESULTS 

In this section we present our main results of con- 
straints on this phenomenological model from the current 
observational data, as shown in Table III. 

In Fig[4] we illustrate the posterior distribution of a 
from the current data. Firstly, we neglect the high red- 
shift probes. The result shows that the current observa- 
tions yield a strong constraint on the parameter: 



a = 0.263 ±0.175 (la) 



(28) 



One can see that the pure ACDM model (a — 0) still 
fits data very well at 2a uncertainty, which is consistent 
with the current status of global fitting results [l|, l32l ]. 
And the 95% upper limit is a < 0.686, which implies 
that there is a significant tension between the flat DGP 
model (a — 1) and the current observations, which is 
consistent with other works (e.g. Ref.[H, Hj). How- 
ever, unlike other works [32], in our analysis we use the 
"WMAP distance information" which includes the "shift 
parameter" R, the "acoustic scale" I a, and the photon 




FIG. 4: One dimensional constraint on the parameter a from 
the different current data combinations: CMB±BAO±SN 
(blue dashed lines) and all real data (red solid lines). 



decoupling epoch z* , instead of R only, to constrain this 
phenomenological model. Recently, many results show 
that the "WMAP distance information" can give the sim- 
ilar constraints, when compared with the results from the 
full CMB power spectrum [35[. By contrast, R could not 
be an accurate substitute for the full CMB data and may 
in principle give some misleading results [3(|. 

And then, we include some high redshift probes, such 
as GRB and LGF data sets. From Table III and FigH 
we can find that the constraint on a becomes slightly 
tighter: 



a = 0.254 ±0.153 (la) 



(29) 



and a < 0.541 at 2a confidence level. As we have men- 
tioned before, the effective equation of state of this phe- 
nomenological model will depart from the cosmological 
constant boundary at high redshifts. Therefore, these 
high redshift observations are helpful to improve the con- 
straints on this phenomenological model. 

From 



These results (Eqs. (|28ll29|) ') are not surprising. 
Fig[T] we find that the effective equation of state of the 
flat DGP model, w c g w — 1 ± a/2 = —0.5, will depart 
from the cosmological constant w = — 1 at high-redshift 
universe significantly. But the current constraint on w is 
closed to w = — 1 0, HH, so we require the small value 
of a to match the current observations. There is a small 
difference that our result slightly favors a non-zero value 
of a, but not significantly, which needs more accurate 
measurements to verify it further. 

In Fig|5]we plot the two dimensional constraint in the 
(£l m ,a) panel. f) TO and a are strongly anti-correlated. 
The reason of this degeneracy is that the constraint 
mainly comes from the luminosity and angular diameter 
distance information. From Eq. ([3]) and Eq. ((4]) we can see 
that when a is increased, the contribution of last a term 
to the expansion rate will become large, due to the pos- 
itive E(z). Consequently, f2 m must be decreased corre- 
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FIG. 5: 68% and 95% constraints in the (Q m ,a) plane from FIG ' 6: 68% and 95% constraints in the (Q m ,a) plane from 



the current observations. 



the future measurements. 



spondingly in order to produce the same expansion rate. 
When combining those current observational data, the 
matter energy density has been constrained very strin- 
gent: Q m = 0.277 ± 0.017 (la), which is also consistent 
with the current status of global fitting results [l], 
Naturally, the constraint on a will also be improved, be- 
cause of the tight constraint on the matter energy density. 

Furthermore, we also investigate the limit on the 
growth index 7 and obtain 7 = 0.570 ± 0.205 at 68% 
confidence level. Obviously, the growth index of the pure 
ACDM 7 = 6/11 0.545 is consistent with this result. 
However, the theoretical value of growth index in the flat 
DGP model, 7 = 11/16 = 0.6875, is disfavored. 

V. FUTURE CONSTRAINTS 

Since the present data clearly do not give very strin- 
gent constraint on the parameter a, it is worthwhile to 
discuss whether future data could determine a conclu- 
sively. For that purpose we have performed an analysis 
and chosen the fiducial model as the mean values of Table 
III obtained from the current constraints. 

The projected satellite SNAP (Supernova / Accelera- 
tion Probe) would be a space based telescope with a one 
square degree field of view with 10 9 pixels. It aims to 
increase the discovery rate for SNIa to about 2000 per 
year in the redshift range 0.2 < z < 1.7. In this paper 
we simulate about 2000 SNIa according to the forecast 
distribution of the SNAP [33|. For the error, we follow 
the Ref.(37j which takes the magnitude dispersion 0.15 
and the systematic error a sys — 0.02 x z/1.7. The whole 
error for each data is given by: 



CTmag(^) = J CT s 2 ys (^) + , (30) 

y 

where n, is the number of supernovae of the i'th redshift 
bin. Furthermore, we add as an external data set a mock 
dataset of 400 GRBs, in the redshift range < z < 6.4 



with an intrinsic dispersion in the distance modulus of 
(7 M = 0.16 and with a redshift distribution very similar 
to that of Figure 1 of Ref.[38|. 

For the linear growth factors data, we simulate the 
mock data from the fiducial model with the error bars 
reduced by a factor of two. This is probably reasonable 
given the larger amounts of galaxy power spectrum and 
lyman-cv forest power spectrum data that will become 
available soon as long with a better control of systematic 
errors in the next generated large scale structure survey. 
In addition we also assume a Gaussian prior on the mat- 
ter energy density fl m as a = 0.007, which is close to 
future Planck constraints [3§| . 

From Table III it is clear that the future measurements 
with higher precision could improve the constraints dra- 
matically. The standard derivation of a is reduced by 
a factor two. Assuming the mean value remains un- 
changed in the future, the non-zero value of a will be 
confirmed around 3cr confidence level by the future mea- 
surements. In addition we also illustrate the two dimen- 
sional contour of parameters Q m and a in Fig(6] Com- 
paring with the contour in Fig|5j the allowed parameter 
region has been shrunk significantly. The future measure- 
ments could have enough ability to distinguish between 
the modified gravity model and the pure ACDM model. 



VI. CONCLUSIONS 

As an alternative approach to generate the late-time 
acceleration of the expansion of our universe, models of 
modifications of gravity have attracted a lot of interests 
in the phenomenological studies recently. In this note we 
investigate an interesting phenomenological model which 
interpolates between the pure ACDM model and the flat 
DGP braneworld model with an additional parameter a. 

Firstly, we find that when a is less than zero, the 
growth function of density perturbation 5(a) will appear 
an apparent singularity. This is because the variable 
will change the sign during the evolution of our universe. 
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And then the (3 term caused by the modified gravity 
model will be divergent at some redshift z t . 

From the current CMB, BAO and SNIa data, we ob- 
tain a tight constraint on the parameter a = 0.263 ± 
0.175 (la), which implies that the flat DGP model (a — 
1) is incompatible with the current observations, while 
the pure ACDM model still fits the data very well. When 
adding the high-redshift GRB and LGF data, the con- 
straint is more stringent a = 0.254 ± 0.153 (ler), which 



means that these high redshift observations are helpful to 
improve the constraints on this phcnomcnological model. 

Finally, we simulate the future measurements with 
higher precisions to limit this phenomenological model. 
And we find that these accurate probes will be helpful to 
improve the constraints on the parameters of the model 
and could distinguish between the modified gravity model 
and the pure ACDM model. 
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